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Non-periodic tilings of Zjv are a difficult to obtain key to a wide range of complex math¬ 
ematical issues. In a musician’s eyes, they are called rhythmic Vuza canons. After a short 
summary of the main properties of rhythmic tiling canons and after stressing the importance 
of Vuza canons, this article presents a new approach: modulo p Vuza canons. Working modulo 
p (especially with p = 2) allows Vuza canons to be computed very quickly. But going back to 
traditional Vuza canons requires understanding of when and how the tiling is a covering. The 
article’s main result is a complete case construction of a modulo 2 tiling. This first solution 
of a modulo 2 tiling is proved in a constructive way, which, it is hoped, is expandable to 
every modulo p tiling. Following this path may lead to a solution to the problem of building 
modulo p tilings and thus a solution to the problem of easily finding traditional Vuza canons. 
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Introduction 


1.1. Context 

The science of rhythmic tiling canons came into being thanks to both Dan Tudor Vuza 
(1991-1993), who approached the subject from a musical point of view, and Gyorgy 


Hajos (1949), who previously developed a mathematical theory -the notion of finite 


abelian group factorization (while he was solving the Minkowski conjecture)- that was 


later linked to Vuza’s work by Moreno Andreatta (1996). 


Musically, a rhythmic tiling canon is a canon where only the onsets matter. A musical 
note is usually described by its pitch, onset, offset, velocity... The canon is said rhythmic 
because most of the note classical attributes are not relevant: it is a canon where the main 
voice is a pattern of onsets. The tiling property of the canon comes from the constraint 
that the pattern of onsets and its translates -the other voices of the canon- must tile the 
time axis: one must hear at every beat one and only one onset. 


Thanks to Vuza’s modelisation of rhythm (Vuza 1985), such a canon is equivalent to 


a tiling of the discrete line Z by a finite integers set and its translates. This, as we will 
explain later, can be represented by a direct sum of a two integers sets being equal to Z. 


Lagarias and Wang have proved in a more general case (Lagarias and Wang 1996 


theorem 5) that a tiling of the line by translations is periodic. This result was already 


known by Hajos and DeBruijn (Hajos 1949), although they never proved it in details. 


* Corresponding author. Email: Helianthe.Caure@ircam.fr 
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Hence, for modelling rhythmic tiling canons, we will not focus on a tiling of Z by a finite 
set of integers, but on a tiling of Zjv- 

For a more complete history of rhythmic canons and recent research on the topic, the 
reader is invited to refer to the special issues on the topic which appeared in both Journal 
of Mathematics &: Music (2009) and Perspectives of New Music (2011). 


1.2. Model 

Definition 1.1 We will denote equivalently a rhythmic pattern by 

(1) a finite set A of non negative integers containing 0 

(2) a polynomial A{X) G {0,1} [X] 

(3) an infinite word with finitely many non-zero letters or a finite word Hq-i over the 
alphabet {0,1}. 

By denoting for any word u, the subword of length k starting at the n-th letter 
u[n,... ,n k — 1] yk G W, the equivalent notations are linked by 

= EaeA X°' and Vn G N Hq-i [n] = where is the characteristic 

function of A and its length verifies maxH + 1 < |^o-i| < oo. 

Elements of a rhythmic pattern are called onsets. 

Definition 1.2 Let ^ be a standard set-theoretic operation (like the sum or cartesian 
product) between sets of integers, and iV G N. We note the operation defined by 
A^]S[B = {a ^b mod N, a G A,b G B}. 

Definition 1.3 Let N G N*, the rhythmic pattern A tiles Z^r if there exists another 
rhythmic pattern B such that one of the following equivalent conditions is satisfied: 

(1) A®n B = Ztv where © is the direct sum (or injective sum ® : A x B ^ Z) 

(2) A{X) ■ B{X) = 1 + X + ... + mod {X^ - 1) 

(3) = Izn where * is the convolution product. 

We say that {A,B) is a rhythmic tiling canon (RTC) ofZ^i, or that the pattern A tiles 
Ztv with the (pattern of) entries B. 

Remark 1.4 It is important to understand that a rhythmic pattern A that tiles Z^r 
doesn’t have to be included in {0,..., N — 1}. Indeed, it is possible that there exists 
a G A such that N < a, because only the elements of the sum A -\- B are projected in 
Zjv. The reader can even prove that from a RTC {A, B) of Ztv, one can built A' verifying 
this latter condition such that [A', R) is a RTC of Z^y. 

Example 1.5 A rhythmic pattern can be written in the following ways 

(1) A = {0,1,4, 5} 

(2) A{X) = l + X + X^ + X^ 

(3) Ao-i = noon or Ao_i = 1100110000 or Aq-i = 1100110... 

It tiles Zi6 with B = {0, 2, 8,10}. 

Remark 1.6 From now on, instead of writing n times a letter o, we will denote aP = 

aa - ■ ■ a. 

An easy way to represent a RTC (A, B), which helps “seeing” the rhythm, is on a grid 
with its x-axis corresponding to the time, and its y-axis the entrance of new voices from 
B playing the pattern A. An onset of a note a + bGA®]\fB of the RTC is depicted by 
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a black square in position (a + 6 mod N, i), with B [i] = b. 

Example 1.7 One can represent the RTC ({0,1,4, 5} , {0, 2, 8 ,10}) of Zig in the following 
way: 



Definition 1.8 A RTC {A, B) is compact if A ©tv S = A © = [0, Ai — Ij exactly, 

without need of projection in'E m- 

Example 1.9 ({0,1,4, 5} , {0, 2}) is compact, whereas ({0,1, 3,6} , {0,4}) is not. 


Definition 1.10 Let k £ W, and A be a rhythmic pattern which tiles Ztv- We note 
A ^ = A ® {0, N, 2N,... ,{k — 1)A1} its k-eoncatenation. 


see: 


As shown by Amiotj (2005) some transformations of RTC give other RTC, 
Proposition 1.11 (Duality) {A,B) is a RTC o/Ztv iff{B,A) is a RTC o/Zjv 


Proposition 1.12 (Concatenation) Let k £ N*, {A,B) is a RTC of "Ln iff {A^,B) is 
a RTC of ZfcTV- 


Example 1.13 The dual of the canon depicted on example |1.7| is the 2-concatenation of 
the dual of the first canon of example |1.9[ 


1.3. Vuza Canons 

Definition 1.14 A RTC {A,B) ofLj^ is periodic if there exists 0 < k < N such that 
A ©TV {k} = A or B ©tv {k} = B. 

If a canon is not periodic, it is called a Vuza canon (VC). 


Theorem 1.15 (Amiot (2005)) Every RTC can he deduced by concatenation and duality 
transformations from VC and the trivial canon ({0} , {0}). 


Example 1.16 {0,1,4, 5} © {0, 2} = Zs is concatenated from {0,1} © {0, 2} = Z 4 , whose 

dual is concatenated from the dual of { 0 , 1 } © { 0 } = Z 2 , concatenated from the trivial 
canon. 


concatenated from ^_I 


concatenated from 



dual of 



concatenated from the trivial canon: 


Remark 1.17 It means that VC are the basis of RTC under concatenation operation. 
Thus, if one manages to obtain every VC, one has a constructive way to get every RTC. 
Note that in practice, it is difficult to built RTC from VC, even if we know that they 
exist and the size of their tile: 
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Theorem 1.18 There exist VC ofLjsi for, and only for N not of the form: 

N = p°,N = p'^q, N = p‘^q^, N = pgr, N = p^qr, N = pqrs 


with p, q, r, s distinct primes. 


Proof. This full result has been shown step by step by Hajos (1949), Redei (1950), 


De Bruijn (1953), Sands (1957) and later by Vuza (1985) independently. 


Those non-periodic canons are also the key to understanding whether a given rhythmic 
pattern tiles or not: 

The tiling condition 


vN _ 1 

A{X) ■ B{X) = 1 + X + ... + X^-^ = mod {X^ - 1) 

A — 1 


gives the idea of splitting the iV-th roots of nnity between the factors A and B, hence, 
looking for their cyclotomic factors. 


Definition 1.19 Let A be a rhythmic pattern, we note 


Ra = {d G N*, the d-th cyclotomic polynomial divides A{X)} , 


Sa = £ Ra,P prime, a G N*} . 

Definition 1.20 Let A he a rhythmic pattern, the following three conditions are called 
Coven-Meyerowitz conditions: 

(To) : A tiles 

(Ti): 71(1) = 0^.65. P 

(T 2 ) : If Pi,P 2 , ... ,p7 G Sa then pf ■ P 2 ' ■ ■ ■ ' Pr G Ra with pi distinct primes. 


They were introduced by Coven and Meyerowitz (1999) in the currently best attempt 
to get a necessary and sufficient condition to know if whether a rhythmic pattern can 
tile or not. 


Theorem 1.21 (Coven and Meyerowitz (1999)) 

(1) (To) ^ (Ti) 

(2) (Ti) A (T 2 ) ^ (To) 

(3) If‘^A has at most two prime factors, then (Tq) =7 (Ti) A (T 2 ) 


Proposition 1.22 (Amiot (2005)) Condition (T2) is closed under duality and concate¬ 
nation transformations. 

Corollary 1.23 We have (Tq) =7 (T2) for every RTC iff (Tq) =7 (T2) for every VC. 

Remark 1.24 This corollary emphasises the importance of VC. Those non-periodic 
canons are nonetheless the construction basis of any RTC, but they are also the bot¬ 
tom line to understand if a given rhythmic canon would tile. 

Though, enumeration of VC is still out of reach of mathematicians: the best we cur¬ 
rently have is an exploration algorithm to obtain slowly all the VC with a given period. 
It was improved recently by Kolountzakis and Matolcsi (2009), but it is still exponential. 
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The next section endeavours to look at VC in a more complex space, which makes them 
easier to obtain, with a still distant hope of going back to VC of Zjv in a non-exponential 
way. 


2. Modulus p tiling 


Many RTC properties are formulated in polynomial notation (like 1.21), and the ring 
Z [V] / [X^ — 1) seems fitted for their study. However, this ring is not factorial, hence, 
one cannot use the decomposition in irreducible elements of the complete polynomial 
1 -|- V + ... -|- to obtain all the factors in {0,1} [X]. 

The polynomial notation of the tiling condition (2nd point or 1.3) gives the idea, 
introduced in ( Amiot|[2004 ) to work in the set {0,1} [X] of polynomials with coefficients 
0 or 1, but it is not a ring. Thus, we try to work in the ring F 2 [X] so that we have 


F 2 [X] X F 2 [X] ^ F 2 [X] 

{A,B) I —> A-B 

and the tiling condition becomes the divisibility of 1 + X + ... -|- X'^"^ by 
A{X) mod {X^ — 1,2), i.e. the classic divisibility condition in Z[X]/(X^ — 1,2) = 
F 2 [X] /(X'^ —1), where (X^ —1, 2) denotes the ideal generated by the set {2, X^ — l} C 
Z[X]. 

We then extend the idea to the ring Fp [X], p prime, and we obtain the notion of tiling 
modulo p. 

Definition 2.1 {A,B) is a RTC o/T^n modulo p if 

A{X) ■ B{X) = 1 + X + ■■■ + X^-i mod (X^ - l,p). 

Modulo p tiling is relevant from a mathematical point of view, but also musically, 
because it enriches RTC with harmony, allowing notes superposition. From now, we will 
denote RTpC for a modulo p RTC. 

Remark 2.2 The equivalent notations with sets from RTC are adaptable with multisets 
to RTpC: 

Let A and B two rhythmic patterns, {A,B) is a RTpC of Z^r if the multiset 
C = A + 7 V B = {a + b mod N,a € A,b € B} verifies lc(^) = 1 mod p if 
n G {0, ...,X — 1}, and lc(^) = 0 elsewhere, with Ic being the multiset character¬ 
istic function. 

Remark 2.3 When we consider a finite multiset C, we can also adapt the notation Cq-i 
to RTpC: 

Let C be a finite multiset of non-negative integers, we note Co-p the infinite word 
with finitely many non-zero letters or the finite word, over the alphabet { 0 ,... ,p — 1 }, 
verifying Vn G N, Co-p [n] = lc(^) mod p and max (7 -|- 1 < ICo-pl < 00 . 

To avoid the confusion between the classical notation Cq-i and Cq -2 for a RT 2 C, we 
do not use the notation C'o_(p_i) for this multiset, even if it seems more fitting. 

Example 2.4 ({0,1, 3} , {0, 2, 3}) is a RT 2 C of Z 7 : 

(1) A+jB = {0,1,2,3,3,3,4,5,61 
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(2) A{X) ■ B{X) = l + X + ...+X^ mod - 1, 2) 

(3) {A + B)o_2 = 1111111 = 111111100000 ... 

Remark 2.5 For any rhythmic pattern A, we have that ^o-i = ^o-p- 

Definition 2.6 Let A and B two rhythmic patterns that could tile or not. Since they 
are non empty (both contains 0), the word {AB)q-p starts with a certain number of 
Is.' the beginning of the tiling. We note m the first index where {A, B) does not tile, i.e. 
the smallest index m such that {A + B)o-p [m] 1 mod p. 

Since both rhythmic patterns are finite, by definition, the word {A + B)o-p is finite with 
(a finite or infinite number of) Os at the end. We note M the last index such that the 
last onset of {A + B) is at time M, i.e. the largest index M such that {A + il)o-p 0. 

Definition 2.7 With previous notations, a modulo p under-cover UCp of {A,B) is a 
finite subsequence of {A + il)o-p of length at least M — m + 1, such that 

UCp [n] = (A + B)Q-p[n + m] 

or equivalently this word starting with a number u of Is, 0 < u < m, or ending with a 
finite number of Os. We say that this under-cover starts at the index i = m — u and when 
relevant, we note UCp{i). 

Example 2.8 Let A = {0,1, 3} and B = {0,1,4}, then A B = {0,1,1, 2, 3,4,4, 5, 7}, 
and a modulo 3 under-cover starting at index 1 is UC^^l) = 2112101 or UC 3 = 
2112101000. 


A + B 



V_^ ^^ 


uc^il) 


Remark 2.9 An under-cover of {A,B) is uniquely defined by its length and starting 
index. 

Its length has to be at least M — m -|- 1 so that the under-cover provides all the 
information needed to fill it (in the sense of definition 2 . 11 ) into an actual tiling. 


_ ^^ 

Remark 2.10 If {A,B) is a RTpC, then any under-cover will be of the form 1 0 , with 
k,l namely the empty word if k = I = 0. 


Definition 2.11 Let A and B two rhythmic patterns, UCp the under-cover of {A,B) 
starting at index i and of length n = \UCp\. 

We say that we fill this under-cover with R' C N finite if 


-I- (R' -I- {i}) ^^ [i,... ,i + n — 1] — 1 . 


Remark 2.12 Note that B' does not depend on B in the sense that if {A, Bfi) and {A, B 2 ) 
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give the same under-cover: UCp{i)i = UCp{i) 2 , we will need the same to fill it. But 
the minimal size of an under-cover depends on both A and B. 


Example 2.13 Let A = {0,1,4} and B = {0,2,5}. We have A + B = 

{0,1, 2, 3,4, 5, 6, 6, 9}, hence {A + B)q _2 = 1111110001 (= 1111112001) and an under¬ 
cover starting at index 6 is [7(72(6) = 0001. We can fill this under-cover with B' = 
{0,2,3}. 

Indeed, we have A + {BU{B' + {6})) = {0,1, 2, 3,4, 5, 6, 6,6, 7,8,9, 9, 9,10,10,12,13}, 
thus {A + {BVJ {B' -I- {6}))) o -2 [6,9] = 1111. 

Note that if we had chosen a longer under-cover, see [7(72(6) = 000100, it would have 
been filled with B' = {0, 2, 3,4} and not B', because {A + {BU {B' + {6})))o-2 [6,11] = 
111111 and {A + {BU {B' + {6})))o-2 [6,11] = 111100 / T®. 


A + B 


A+(ba[ +{ 6 })) 

A + {B U {B'+ {6})) 



A + B' 




0 0 0 +1 
A + W 




0 0 0+10 0 


Remark 2.14 When one wants to tile with a pattern A, one can imagine filling the null 

—k 

under-cover 0 for some k G N, starting at index 0. 

This is the meaning of an under-cover: it can be seen as a partial covering that one 
needs to fill to obtain a tiling. This is how the greedy algorithm presented below works. 


Definition 2.15 If one fills an under-cover UCp{i) of {A,B) with B', then we call 
‘following under-cover” an under-cover of {A, B U {B' + {*})), starting at index i + 
\UCp{i)\. It will be denoted UCp. 


Example 2.16 In 
17(7^(10) = 0011. 


the previous example 2.13 


a following under-cover of [7(72(6) is 


Definition 2.17 {A,B) is a compact RTpC of'E m if 


A{X) ■ B{X) = l + X + ... + X^-^ mod p. 


Theorem 2.18 (Warusfel 1971) Let p prime and P G¥p [77] such that P{0) + 0, then 
there exists N G N such that P divides X^ — 1. 


Remark 2.19 As hrst noticed by Amiot (2005), this theorem entails that any rhythmic 
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pattern tiles modulo p: 

Let A{X) G Z [X] be a rhythmic pattern and a fixed p prime. Considering the polyno¬ 
mial PiX) = A{X){X — 1) G Fp [X], and the quotient polynomial Q G Fp [X] given by 
the theorem such that P{X) ■ Q{X) = (X^ — 1) mod p we can obtain the tiling condition 
P(X) • Q{X) = [X^ — 1) mod (X-^ — l,p) with Q a polynomial representing a rhythmic 
pattern. Indeed, by iterating the transformation 

aX^ = {a- l)X^ + X^+^ mod (X^ - 1 ) 

over the polynomial Q, it can be associated to a rhythmic pattern Q G Z[X], i.e. a 
polynomial of {0,1} [X]. 

Note that in F 2 [X], every polynomial has its coefficients in {0,1} which means that 
this transformation is unnecessary, hence every rhythmic pattern gives a compact RT 2 C. 
Indeed, it is clear that the non-compacity of a RTC (in Z or Fp) comes from the projection 
modulo (X'^ — 1); i.e. here from the transformation we apply on Q. 

Example 2.20 In F 3 [X], let ^(X) = 1 -|- X -|- X^ be a rhythmic pattern. It verihes 

x^^ - 1 = ^(X)(X - 1) • (X2° -h X^® -h 2X^^ -h 2X^® -h 2X^^ + 2X^^ 

+2X^^ -h X^ -h 2X® + X'^ + X^ + X^ + 2X^ + 2X3 x‘^ + 1) 

The quotient polynomial Q(X) = X^o + X^^ + 2X^’^ + 2X^^ + 2X^^ + 2X^3 2 X ^0 + 

X® -|- 2X3 _|_ _j_ -j- + 2X‘^ -|- 2X3 _|_ _|_ 2 (ioes not have its coefficients in {0,1}, 

thus is not a rhythmic pattern. By applying the transformation 2X^ = X^ + j^fc+24 
its terms (2 * X^\ 2 * X^^ 2 * X^^ 2 * X33, 2 * X^o, 2 * X®, 2 * X^, 2 * X3), we obtain the 
polynomial 


Q{X) = X^O + x39 + X^3 + x^'^ + X39 + x33 + X33 + X^^ + X37 

+X^3 ^ ;sj:34 ^ + X32 + X® + X^ + X® + X® + X^® 

+X^ + x27 + x3 + X V 1 


which is a rhythmic pattern verifying 

A{X) • Q(X) = 1 + X + ... + X23 (j^24 _ 3 ^^ 


This theorem has inspired a modulo p tiling algorithm in (Amiot 2011). We present 


here an algorithm devised from the previous one, but improved inasmuch as it is greedy 
and optimal for any rhythmic pattern which gives a compact RTpC, a fortiori for any 


rhythmic pattern with p = 2 thanks to remark 2.19 
Require: A C N finite and 0 G A 


B = {0}, N = max A 

while (A -|- B)o-p 7 ^ 1 for some n do 

i <— first index such that (A -I- R)o-p 7^ 1 mod p 
if i 7 ^ (maxR -I- max A -|- 1) then 
B := BU{i} 

N := i + max A 
else 

break while 
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9: end if 

10 : end while 
11 : return {B, N) 

Example 2.21 This algorithm is really easy to understand from a graphical point of 
view. Indeed, if one considers the representation in a grid, it means adding the rhythmic 
pattern A at the earliest time when there is not 1 mod p onset. It is also a good way to 
understand the notions of definitions 12.Ill and 12.151 
See the algorithm applied to = {0,1,4}. For better understanding, at every step 
depicted of the algorithm, we will show the value of [A + i?)o -2 and of the smaller 
uner-cover [7(72 (i) the algorithm is trying to fill. 


{A + B)o_2 = 11001000 -•• 

17(72 ( 2 ) = 001 


(71 + .B)o-2 = 1111101 
17 ( 72 ( 5 ) = 01 


{A + B)o-2 = 1111110001 
17(72(6) = 0001 


{A + B)o-2 = 11111111011 

17(72(8) = on 


Step 1: 




With i = 6: 

















1” 1 

l" 


■1 


1 ' 




With i = 10: 
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{A + B)o-2 = 111111111111111 

1/(72(15) = empty word 


Remark 2.22 Such a tiling is not unique: the rhythmic pattern A = {0,1,4} tiles with 
B = {0, 2, 5,6, 8, 9,10}, but also with the following Bi’. 

• Bi = {0, 2, 5,6,8, 9,10,15,17, 20, 21, 23, 24, 25} 

• B2 = {0, 2, 5,6,8, 9,10,15,17, 20, 21, 23, 24, 25,30, 32, 35,36, 38, 39,40} 


Proposition 2.23 This greedy algorithm is optimal in the sense that it gives the small¬ 
est number of entries ^B (and so the smallest N) such that (A,B) is a compact RTpC 
o/Zat. 

Proof. Let p be a prime number, and A a rhythmic pattern such that it admits a compact 
RTpC. Reductio ad absurdum, let us suppose there exists B such that {A, B) is a compact 
RTpC of Zjy and jJR < (jR with B the output of the algorithm. 

B B because since {A, B) is a compact RTpC, it means that during the algorithm, 
the smallest element i oi B \ B would verifie line 4: i = maxR + max^ + 1 and hence 
B = B. 

We can suppose without loss of generality that B and B are sorted. Let i the smallest 
index such that B[i] = b < B[i] = b. 

After i — 1 loops, the algorithm produces the subset R [0,... i — 1] providing the same 
under-cover U (7p(0) as the subset R [0,... i — 1]: 

UCp{0) = (A+R [0,..., i — l])o-p [0,..., iV — 1] = (A+R [0,..., i — l])o-p [0,..., — 1]. 


Since R tiles with A, at time b, there cannot be exactly 1 mod p onset in the under¬ 
cover: UCp{0) 5/1 mod p. Step 3 of the algorithm, during the i-th loop, returns the 

smallest index V such that at this index of the under-cover there is anything but 1 mod 
p onset, and will give a b' < b < b, which is a contradiction with the fact that the Lth 
element of R is b. 


Example 2.24 As we can see in example |2.21 
modulo 2 with A = {0,1,4} is 7. 


the smallest number of entries to tile 


Remark 2.25 Thanks to theorem 1.15, we have seen that VC are the basis of RTC under 
the concatenation operation. Also, corollary 1.23 gives us the result that VC are the key 
to attaining a necessary and sufficient condition on the tiling property of a rhythmic 
pattern. 

This greedy algorithm provides minimal RTpC under concatenation. They are the 
“modulo p Vuza canons”. We hope that a connection can be established between both 
kinds of Vuza canons, the latter modulo p being relatively easy to compute, to be able 
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to return to classical VC and thus obtain a faster way to enumerate them. It remains, 
however, a distant ambition. There are two principal differences between classical VC and 
the RTpC ones we are working with. First, the algorithm only computes compact RTpC, 
whereas VC are not usually compact. Future research, with the aim of linking back those 
two types of canons, will have to search and compute minimal non compact mod p Vuza 
canons, i.e. minimal under concatenation and duality, but also tiling a smaller period 
than the compact one. The second difference is of course the distinction between tiling 
Z and tiling Fp. A VC admits one and only one onset every beat, unlike RTpCs which 
admit one modulo p per beat. So we need to understand how to remove those extra 
onsets. The paper will now focus on a better understanding of this latter challenge. 

The goal is to understand when and how we are able to obtain ‘donsets’, a neologism 
introduced here to express the fact that the onsets have copies at this time, giving us a 
covering instead of a tiling. This term is also a contraction of d-onset, and if we have 3 
layered notes, we can talk about a 3-onset. 

Definition 2.26 We note the multiset of donsets of the RTpC {A,B) ofL^: 

D(a,b),p = IJ U 

nGZiv kp^ 

fcGN", 

l(A + „i3)W = fcP+l 


One can extend the notion of duality to compact RTpC: 

Proposition 2.27 {A,B) is a (resp. compact) RTpC iff {B,A) is a (resp. compact) 

RTpC and 

Proof. Let (A, B) be a (resp. compact) RTpC of Z^r. C = A Tw B (resp. C = A + B) is 
a multiset such that: 


C = 


0,0,...,0,1,1,...,1,...,A^-1,A^-1,...,V-1 


kop-\-l times fcip+1 times 


/cat-iP+I times 


with Vi G [0, V — IJ, ki £ N. 

By commutativity in Z, we get that C = B +]\f A (resp. C = B + A), hence {B, A) is 
a (resp. compact) RTC of Zw modulo p, and the donsets correspond to the ki 0 and 
are at the same positions a + b mod N = b + a mod N (resp. a + b = b a) and in the 
same amount kip. ■ 

Remark 2.28 However, if B is the smallest rhythmic pattern obtained by the algorithm 
for a given A, A may not be the smallest complement obtained by the algorithm if the 
entry is R, if A is a concatenation of different rhythmic patterns that tile with B. 

For instance, if the rhythmic pattern in entry is A = {0,1, 2, 6, 7, 8}, the algorithm 
modulo 2 returns B = {0, 3} for which it returns A' = {0,1,2}. In such cases, we always 
have that such an A' is always a subset of A. 


Remark 2.29 If the tiling is not compact, the best way we currently have to compute 
the size of a RTC quickly is the algorithm presented in (Amiot 2011) which returns a 
multiple of the size of the smallest RTC with a given rhythmic pattern; this is achieved 
by incrementing the suitable values of N until we can find one such that A{X) divides 
— 1 modulo p. One must bear in mind that in order to discover the exact size of the 


11 










August 25, 2015 Journal of Mathematics and Music VuzapHCreviewS 


n 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

N 

3 

7 

15 

21 

63 

127 

63 

73 

889 

1533 

3255 

7905 

11811 

32767 

255 


Table 1.; Size N of the minimal compact RT 2 C with the rhythmic pattern A = {0,1, n}. 


smallest RTC, one needs the greedy algorithm when the tiling is compact. 

See, in Table the minimal size of RT 2 C (hence compact) of rhythmic pattern A = 
{ 0 , l,n} for some small n obtained with the greedy algorithm. 

One can notice that rhythmic patterns of the form Ak = { 0 , 1 , 2 ^| give smaller compact 
RTC. We will now focus on these patterns. 

3. The case = |0,1, 2^} 

The following theorem is the main new result of this paper, and is the first cardinality 
result for RTpC. 

Theorem 3.1 For all k G N*, the rhythmic patterns couple {A^, B^) = (|0,1,2^} , Bj.) 
is the smallest compact RT 2 C of Zjv, and is of size N = 4^ — 1 with ^B^ = 4^ — 3^ 

3fc+i _ 

entries, and ttR*(Afc,Rfc ),2 = 4^-^- donsets. 

Thanks to a constructive proof, this theorem allows us to build the pattern of entries 
Bk which tiles modulo 2 with a given A^ = { 0 , 1 , 2 ^|. Since one can also easily prove 
that in this case one cannot have d-onsets with d > 3, one arrives at a comprehensive 
enumeration: the number of entries, the size of the RTC, and the number of donsets. 

Before giving a formal proof, we will start with an explanation of how those tilings 
are built, which will in turn demonstrate the general concept of the proof. Some tech¬ 
nical lemmas will form the second part of this proof, and we will hnish with the formal 
demonstration of the theorem. 


3.1. Introduction 


Let us begin with a close look at how the greedy algorithm builds a tiling with the 
pattern A^ = { 0 , 1 , 8 }. 

At step 1, the algorithm places the pattern A 3 at the first position, that is to say 
B = {0}, hence (A -h R)o -2 = 110000001. 

The algorithm then enters the while loop, and tries first to fill the under-cover 
C/C' 2 ( 2 ) = 0000001. R is filled with the entries B' = {0,2,4} which gives after three 
steps 3 of the algorithm B = B U ({0,2,4} -|- {2}) = {0,2,4, 6 }, thus (A -|- R)o -2 = 
111111111010101. The algorithm goes again in the while loop, trying now to hll the 
under-cover 17 ( 72 ( 9 ) = 010101. And so on. 

The idea is to fill the current under-cover, and doing so, we get a new under-cover: 
the “following under-cover” of definition 2.15 ; and we keep on hlling it until we obtain 


a compact tiling, that is to say an empty or null under-cover. The entries pattern is 
built in this incremental way, filling the following under-covers one after another during 
the while loop. 


Remark 3.2 Let i G N, our rhythmic pattern A^ = {0, 1 , 2 ^} being of size max{Ak) — 
min{A}f} = 2^, if, to fill any under-cover of size 2^ starting at index i (and ending at 
index z -|- 2^ — 1) we need the entries B', then the following under-cover of size 2^ (and 
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you can find one of this size) and starting at index i + 2^ will be Bq_i. This fact is due 
to the gap between the first onsets {0,1} and the third and last onset {2^} of ^4^. 

Example 3.3 Let A 2 = {0,1,4} and B = {0,2} which fills the null under-cover 0000 
starting at index 0. We then obtain (A 2 + B)q -2 = 11111010, hence the following under¬ 
cover starting at index 2^ -|- 0 = 4 and of size 4 is 1010 = Bq-i. 


A + B, Bo_i = 1010 


A + B', B'q_^ = 0110 


A + B", Bi;_i = 1110 


0000 





An efficient way to exploit this phenomenon is to built Bj. in a recursive way. If we 
write (i?fc)o-i in an array V{k) with 2^ columns, from left to right then top to bottom, we 
obtain that every line I oiV{k) is the following under-cover of the under-cover obtained 
with Afc and (i?fc)o-i read in lines 1... ^ — 1. 

Example 3.4 For a better understanding, here are A^, Bj., (ilfc)o-i and (ilfc)o-i in the 
V{k) array for /c = 1, 2, 3: 

. k = l,Ai = {0,1, 2},Bi = {0} , (5i)o-i = 1 

F(l) = (10) 

. k = 2,A2 = {0,1,4} , ^2 = {0, 2, 5, 6, 8, 9,10} , 

(52)0-1 = 10100110111 
/l 0 1 0\ 

F(2) =0110 

Vino; 

• A: = 3, As = {0,1,8}, 

53 = {0, 2,4, 6, 9,10,13,14,16,17,18, 20, 21, 22, 27, 28, 29, 30,32, 

34,35,36,37,38,41,42,43,44,45,46,48,49,50,51,52, 53,54} 


(. 83 ) 0-1 


F(3) 


= 10101010011001101110111000011110101111100111111011111110 
/l 0 1 0 1 0 1 o\ 

01100110 
11101110 
= 00011110 
10111110 
01111110 
Vl 1 111 1 1 0/ 


Definition 3.5 Let k G N*, B a pattern, and UC 2 an under-cover of {Ak,B) of size 

h k 

n <2^. We denote —)• a postfix function over and to the finite words over the alphabet 

k 

{0,1} defined by UC 2 —^ (-B')o-i iff the pattern B of size n (you can pad in some O 5 
in the end to adjust the size) fills UC 2 with A^. 


Example 3.6 


We have O^” 


lO*" if 2n < 2*^. 


See in example 3.4 above, for k = 3 and n = 4, 


one can see that the following under- 
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_g _^ 

cover of 0 of size 8 (which is the first under-cover of the tiling) is 10 . 


Remark 3.7 Thanks to remark 3.2, we know that the pattern such that {Ak,Bk) is 


a RT 2 C can be written in the form {Bk)o-i = Bi ■ ■ ■ Bn with 


\/i < n IBA = 2^ and 0^ 


Bi 


Bn 


and the theorem’s proof will yield n = 2^ — 1. It means that ever y und er-cover Bi will 
be the following under-cover of the under-cover Bi-i as defined in 2.15 


Written in the array V{k) as defined earlier, it means that V{k) = 


(Bi\ 

B 2 

\BnJ 


Example 3.8 With k = 2, we have 0^ 1010 0110 1110 and indeed we have 

( 52 ) 0-1 = 10100110111 = 101001101110 that is to say: 

/ 1010 \ 

V{2) = 0110 . 

V 1110 / 


Definition 3.9 Let us define the array T{k) by induction over k G N*; 

Let T(l) = (10). If one has the array T{k) of size {2^ — 1, 2^) , one can build the array 
T{k + 1) of size (2^+^ — 1,2*’+^) in this way: 


(T{k) 

T{k) 



11... 

11... 

1 0 

1 0 

\ 

11... 10/ 


— u h _ (2^ _li 

with T{k) being an array of size (2^,2^) defined by: its first line is 0^ 1, and its 

last lines are T{k) whose last column o/O is changed into a column of 1 (one can easily 
prove by construction that the last column ofT{k) is null for all k). 


Example 3.10 See for some values of k the array T{k) with separation lines from the 
inductive construction: 



r(3) = 



0 

1 

0 

\ 




0 

1 

1 

0 





1 

1 

1 

0 

) 





0 

1 

0 

1 

0 

1 

0\ 

0 

1 

1 

0 

0 

1 

1 

0 

1 

1 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

1 

0 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

0/ 


where we can see T(2) = 


/O 0 0 1\ 
1 0 1 i 
0 11 i 

villi/ 
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Remark 3.11 The proof of the theorem will be to establish the equality V{k) = T{h) 
for all k € N*. For this purpose we need some technical lemmas presented in the next 
part of this proof. 


3.2. Lemmas 

Lemma 3.12 For all n G N*, to fill an under-cover UC 2 {i)n = Ol^O starting at index i 
with the pattern 11, one needs the entries = 1 0. 


Proof. By induction over n. 

If n = 1, to fill the under-cover UC 2 {i)i = 010 with the pattern 11, one needs to put 
an onset at index i, hence = 1... The new under-cover is 100 = 1(1 -|- 1 mod 2)0. 
One then needs to put the pattern at index i -|- 1 to fill the new 0, hence Bq_^ = 11.. 
The new under-cover being 111, the under-cover is a tiling and we do not need any more 
entry, hence = 110. 

Assume the lemma is true for some n, let us see how one can fill the under-cover 
UC2{i)n+i =01 0 starting at index i with the pattern 11. 

The under-cover starts with a 0, hence the word of the filling pattern starts with a 1, i.e. 
Bq_.^ = IT”'"'"^, and the new under-cover is 10T”0 = 1(1 -|- 1 mod 2)l”0 = 1 • UC 2 {i + l)n- 
By induction hypothesis, to tile the n -|- 2 last index of this under-cover, the n 2 last 
terms of the entries have to be l”^^0. Hence we have Bq_i = ll"’"*~^0 = l”^^0. ■ 


Lemma 3.13 




0(2-1) it(2^-i)0 


Proof. Let us denote 1^^ ^^0 1^^ ^^0 C. 

The first block of Is of size 2^ — 1 is already a tile, there is no need of entries for those 
2*^ — 1 first indexes, hence the word C starts with 0^ 

We then arrive to the second block 0 1^ T, whose length verifies |0 1' T| < 2^~^^. 
Only the first two onsets of = |0,1,2^+^} are used to tile the under-cover, hence 
we try to tile the rest of C with the pattern {0, l}o_i = H- So the filling verifies the 


conditions of lemma 


3.12 


meaning C \2f ... = 1^ 0. 


We indeed have C ^Hl^^ ^^0. 


Lemma 3.14 V/t G N*,Vn < 2 ^ - 2,0l”0 


Proof. Since |0T™0| < 2^, we only tile with the first two onsets 11 of the pattern A^, and 


we are again in the conditions of lemma 3.12 


Lemma 3.15 \/k G N, 01^^ ^^0 


^(20 


3.12 


3.12 


_ (nk _o') K 

Proof. In this case, we don’t have |01’ T| < 2^ and we cannot use lemma 

However, we have exactly |01^ "^01 = 2^. It means that we are out of lemma 

conditions only for the last index. Hence if we have to add a pattern at the first index 

(and since the under-cover 01^ T starts with a 0 we will have to), we can consider the 

—( 2 ^— 2 ) 

same under-cover with an additional onset at its last index 01 1 that we tile with 

the pattern 11. 

By an induction like the one used to prove lemma 


3.12 


—( 2 ^ — 1 ) 

we obtain the entries 1 0. 
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But the 0 on the last index is already filled by the last onset of A^, so we only need to 
add a last time the pattern at the last index, and the following under-cover is 1' . 


Lemma 3.16 VA: G N, if |yl| = \B\ =2^ — 1, if A - 0 —^ B • 0 then A • 1 —^ B • 1, with 
U • u being the concatenation of the word U to the letter u. 

Proof. By definition of an under-cover, adding some zeros at the end does not change it, 
so A = yl • 0. If we denote B (resp. Bl) the rhythmic pattern such that A = A-t) (resp. 
yl • 1) is the UC 2 {i) of {A^, B) (resp. of {A^, Bl)), and associate the words A, B with the 
sets A, B we have by dehnition 


(^Ak + (^B U (B -|- 2 


0,...,f + 2^-l 


= and {Ak + B)o -2 = 1" -A-O. 


Denoting by C* a multiset of the form 


C=( 


0,0,... ,0,1,1,..., 1,... ,f - l,f - 1,... - 1 


2/co+l times 2Aii+l times 


2fci_i+l times 


with Vn G [0, f — Ij, /c^ G N such that Ak + B = C U {A + {i}), we have 


{[C U (A -I- {i})] U [Ak -I- (B -I- {i})])o_2 


0,..., i + 2^= - 1 




or, since \B\ =2^ — 1, 

(C U (A -I- {i}) U ({0,1} -I- (B + {f})))o_2 


0,..., f + 2^^ - 1 


^ Y*j|A-0| 


—2—I A-Ol 

by removing the i-th first indexes of 1 1 that come from C 


((A -I- {i}) U ({0,1} -I- (B -I- {f})))o_2 


0,..., f + 2^= - 1 




SO, by excising the i-th first beats: 


(A U ({0,1} -I- B))q_2 


0,...,2'^-l 


= 


When one fills the under-cover A ■ 0, one needs the entries B ■ 0. Likewise, since the sets 
associated with the words A - D, B ■ 0 are the same as those associated with the words 
A, B by definition, we have, using the same notation 


(A • 0 U ({0,1} -I- B • 0))g_2 



j|A|+l 


The sets A • 1 , B • 1 associated with the words A ■ 1, B ■ 1 verify A-1 = Au{2^}, 
B • 1 = BU {2^}. 

We have 
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(A • 1 U ({0,1} + B • 1)) = ((A U { 2 "}) U ({0,1} + (B U { 2 "}))) 

= ({ 2 '=} U A U ({0, 1} + B) U ({2^ 2^ + 1 }))) 

= ((Au({0,l} + B))U({2^2^2" + l}))) 

Thus 

(A • 1 U ({0,1} + B • l))o_2 = + 1 + 1)1 = 

So, by adding the first i beats with a multiset Cl such that A^+Bl = C'1U(A • 1 + {f}), 
we have 


([Cl U (A • 1 + {f})] U [Ak + (B • 1 + {i})])o_2 


0,...,i + 2^-l 




which means 


Ak+ [BIU {B ■ 1 + {i 


0-2 L 


0,..., f + 2^^ - 1 




we indeed have that A ■ 1 —;■ B ■ 1. 


Lemma 3.17 Vn,/c G N,Va G {0,1}, if |A| = \B\ = 2^ — 1, if A ■ a —^ B ■ 0 then 

-i- n k+n—l ^ 

A • a —)• B ■ 0 . 


Proof. The fact that B • 0 ends with a 0 means that to fill vl • a = UC 2 {i), one does not 
need to put the pattern 11 (first two onsets of A^) at index i + 2^. Hence we do not brim 
over the 2^-sized following under-cover when filling A • a with B. To fill the following 
under-cover we do not need to keep the memory of the previous filling (except for the 
last onset of Ak which creates the following under-cover of course). We can concatenate 
the filling pattern if we want to fill the n-concatenation of the under-cover, taking care 
to tile with |0,1, instead of {0,1,2^} to make sure that we will not run into 

the last onset of A^+n-i and that we will only use the first two onsets 11 during the 
filling. ■ 

Those previous lemmas help for the recursive proof of the equality V{k) = T{k), 
whereas the next one will be used for the cardinality result. 

Lemma 3.18 If one constructs a minimal RT 2 C with A such that ^A = n, then every 
k-onset verifies k < n. 

Proof. Let H = (oi,..., a„}, H such that {A, B) is a RT2C of Zat and suppose a /c-onset 
at time t with k > n. 

This A:-onset is built from different layered onsets of A: ,.. .aj^, and we note bj. the 

voice entry in B that produces the onset aj. played at time t. 

Since /c > n, by the pigeonhole principle, there exists j* / ji such that aj. = aj,. Hence 
the rhythmic pattern is played twice at the same entry bj. = = t — . So, at every 


17 


August 25, 2015 Journal of Mathematics and Music VuzapHCreviewS 


time bj. + am = + Qm, m = 1,..., k, there is two onsets layered from A, that counts 

as 0 mod 2. 

So, if we will have the exact same number of onsets mod 2 at every time oi A B 
and A+n {B\{bj.,bj^}). 

Hence {A, B \ {bj. , bj^}) is a smaller RT 2 C of Z^r, which is a contradiction with the fact 
that {A,B) is minimal. ■ 


3.3. Proof of the Theorem 


To prove theorem 3.1, we will start to prove the 


Theorem 3.19 The array T{k) built previously in definition 3.9 is equal to the array 
B2 


V{k) = 


\BnJ 


7t2^ 


with = Bi • • • Bn, i.e. 0 


k j-y k 

- > Bl -; 


Bn and IBA =2^. 


Proof. Let us prove this equality by induction over k: 

If A; = 1, we have indeed r(l) = (10) = H(l) and the pattern of entries is Bi = {0}. 

If the equality is true for some k, let us prove it for A: + 1: 


By lemma 3.17 with n = 2 we obtain that if T{k) gives the pattern of entries to tile 


with the pattern Ak, then (T(A:)|T(A:)) gives the begining of the pattern of entries to tile 
with Ak+i. Indeed, we proceed by blocks of size before we come to the last onset 
of Ak+i. 


By lemma 3.13 we obtain the construction of the first line of the bottom block of 


T{k + 1), i.e. the first line of T{k) which is 0^^ adjoined to 1^ ^0. 


By lemma 3.14, we obtain the block at the bottom-right of T{k + 1): 

11...10 

11...10 


11...10 

Indeed, one can notice that we always have a 0 before those lines of Is, so we are in the 


condition of lemma 3.14 if we didn’t, it means that we have a line full of Is, and hence 


we are in the case of lemma 3.13 which entails having a complete tiling with Ak+i and 


finishing the equality, see later when we use lemma 3.15 (or going towards step A: -|- 2) 


Lemma 


3.16 


explains why all but the first lines of T{k) look like T(k). We try to tile 
with the pattern 11, since the last onset of A^+i is not taken into account in a block of 
such size. With the recursion hypothesis V{k) = T{k), the procedure gave us the word 
in T{k), and we can think of it as though we were filling a previous line of Os, as it is 

highlighted in remark 


2.14 


_ 2 ^_ 1 

Here the previous line is 0 1, the first line of T{k), hence 

the following word obtained is T{k) but with a 1 instead of a 0 in its last column as 
shown in lemma [3. 161 
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Finally, lemma 


3.15 


—— 1 ) 

proves the end of construction. The last line of T{k + \) is 1^ 


which gives a tiling modulo 2 with Ak+i = {0,1,2^+^} because we use in last position 
(index 2^+^ of this line) the last onset of Indeed this sequence of 2^+^ —1 consecutive 
entries gives a tiling modulo 2 which is followed by — 1 onsets not donsets, finishing 
the tiling. 


We can now conclude by induction that Vfe G N* the array T{k) built like in definition 


3.9 is equal to V{k), i.e. is the array containing (ilfc)o-i written line by line. 


Now that we have proved this hrst result, we know how to construct the pattern of 
entries and in order to obtain the cardinality results in theorem 3.1, we just have to 
tally in T{k). 


Proof. The number of entries b{k) = is the number of Is in the array T{k) and we 
can count them by induction: 

We have 6(1) = 1;6(2) = 7. 

By construction of r(/c + 1), we have the equality 

6(A; + 1) = 3 X 6(A;) + 4^ 

Indeed, the term 3 x h{k) comes from the array T{k) put three times in T{k + 1) and 
the term 4^ is from the bottom right block of Is of size (2^, 2^). 

One can verify easily that for all /c, b{k + 2) — 7b{k + 1) + 12b{k) = 0, hence 

4 3 


To obtain the size N of the RTC, we know that it will end with putting the pattern 
Ak on the last 1 of (Bfc)o-i, so 

iV = (2*^ - 1) X 2^= - 1 + 2^= = 4^^ - 1. 


The hrst product comes from the size (2^, 2^ — 1) of the array T{k), and the next —1 
from the fact that the last 1 in T{k) is on the penultimate position. The last +2^ is from 
the addition of the pattern Ak at this index. 


Finally, to enumerate the donsets, we simply use lemma 3.18 In our case, we have 
‘^Ak = 3, so donsets cannot have more than 3 onsets layered. And since we obtain a 
RT2C, it means that there is 1 or 3 onsets at each time, so every donset is a 3-onset. 
Hence, if we note Dk the set of index where there is a donset, we have the formula 


ttAfc X tIRfc - 2’^Dk = N. 


Indeed the factor 2 multiplying comes from this lemma 3.18, saying that every 
donset is exactly two onsets above a classical tiling; the result follows: jJH/j = — 

3fc+i _ I 
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4. Conclusion 


Finding a way to obtain all the VC of a given period Z^v is an issue for composers. It is 
also a major challenge for mathematicians, as solving the Coven-Meyerowitz conjectnre 
could help to solve the spectral conjecture in dimension 1. This conjecture, also called 
the Fuglede conjecture, has recently been proved to be false for dimensions higher or 
equal to 5 (Tao 2003), and from this work, Matolcsi (2005) has proved it to be false in 
dimensions 3 and 4. The problem remains open in smaller dimensions. 

During the last decade, many approaches have been tested to built and understand 
VC: ex haustively (jKolountzakis and Matolcsi 2009), by using perturbations of periodic 
tilings ( Kolountzakis et al.||2003 ), through transformations that conserve the property of 
being a VC ( Amiot||2005 )... 

This article introduces a new approach with the modulo p VC, and especially modulo 
2, which are substantially easier to compute. It is also a sub-problem of a cryptographic 
issue: factorisation of polynomials in F 2 [V], which is equivalent to the factorisation 
of some multisets. The interested reader can find a recent work on factorisation of a 
multiset into two multisets (Neumann and Praeger 2004). The modulo p VC problem 
can be rephrased as the still unresolved mathematical topic of factorisation of a multiset 
into two proper sets: a new mathematical branch to explore, rising from a musical bough. 
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